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We report and discuss case study simulations of the Rayleigh- Taylor instability in the Boussinesq, 
incompressible regime developed to turbulence. Our main focus is on a statistical analysis of density 
and velocity fluctuations inside of the already developed and growing in size mixing zone. Novel 
observations reported in the article concern self-similarity of the velocity and density fluctuations 
spectra inside of the mixing zone snapshot, independence of the spectra of the horizontal slice level, 
and universality showing itself in a virtual independence of the internal structure of the mixing zone, 
measured in the re-scaled spatial units, of the initial interface perturbations. 
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I. INTRODUCTION 

The Rayleigh- Taylor instability occurs when a heavy 
fluid is being pushed by a light fluid. Two plane-parallel 
layers of fluid, colder on top, are in equilibrium while the 
slightest perturbation leads to the denser fluid moving 
down and the lighter material being displaced upwards. 
The early, linear stage of the instability was described by 
Rayleigh [l[ and Taylor (2j, and summarized in ||. Fur- 
ther development of the instability leads to enhancement 
of the mixing and to a gradual increase of the mixing 
zone, which is the domain where proportions of heavy 
in light and light in heavy are comparable. Dimensional 
arguments, supported by large-scale modeling 0,0], sug- 
gest that the half-width of the mixing zone, h, grows 
quadratically at late time, h oc aAgt 2 , where A is the At- 
wood number characterizing the initial density contrast, 
g is the gravitational acceleration, and a is a dimension- 
less coefficient. 

The coefficient a was the focus of almost every pa- 
per written on the subject of Rayleigh- Taylor turbulence 
(RTT) during the last fifty years. The first attempts to 
look inside the mixing zone were initiated only in late 
1990s @, Uli due to advances in experimental and nu- 
merical techniques. The results of many studies and the 
controversies surrounding the a-coefficient were recently 
summarized in the review combining and analyzing the 
majority of existing a-testing simulations and experi- 
ments @- In this article we also discuss the developed 
regime of RT turbulence. Our main focus is on the anal- 
ysis of the internal structure of the mixing zone, and we 
trace the ev-coefficient only for validation purposes. 

Our analysis of the mixing zone develops and ex- 
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tends previous exp erimental [6|, [lfj, [l2( and numerical 
[E EE EE EE EE EE El observations on the subject, 
and it is also guided by phenomenological considerations 
discussed in [HI]. The essence of the phenomenology, 
which utilizes the classical Kolmogorov-41 approach [2fJ |. 
can be summarized in the following statements: (i) The 
mixing zone width, h, and the energy containing scale, 
i?o, are well separated from the viscous, r\, and diffu- 
sive, rd, scales. In the inertial range, realized within 
the asymptotically large range bounded by Ro/rj from 
above/below, turbulence is adjusted adiabatically to the 
large-scale buoyancy-controlled dynamics, (ii) In three 
dimensions, the velocity fluctuations at smaller scales 
are asymptotically decoupled from weaker buoyancy ef- 
fects (29[. (iii) Typical values of velocity and density 
fluctuations scale the same way as in the stationary, ho- 
mogeneous Kolmogorov turbulence, Sv r ~ (er) 1 / 3 , and 

Sp r ~ ej' /2 e _1 / 6 r 1 / 3 , where the energy Kolmogorov flux, 
e, increases with time while the density fluctuations flux, 
e p , remains constant, according to the buoyancy pre- 
scribed large scale dynamics. All of these three theses of 
the phenomenology are consistent with available experi- 
mental [E E3] and numerical [IE EE EE observations 
of the velocity and density spectra. One particularly im- 
portant consequence of the phenomenology, the decrease 
of the viscous and dissipative scales, was also predicted 
in [l5[ and numerically confirmed in [IE, EU . 

In spite of its relative success in explaining RTT, the 
phenomenology [IE| is, obviously, not free from deficien- 
cies. First, the asymptotic, large time character of the 
theory turns into a handicap in explaining numerical and 
experimental data, taken at finite, and actually mod- 
est, times. Second, the phenomenology treats all z-slices 
within the mixing zone equally. Third, the phenomenol- 
ogy does not differentiate between the mixing zone width, 
h, and the energy containing scale, Rq, for the turbulent 
fluctuations. 

Improving the phenomenology from within itself, or 
by some complementary theoretical means, does not 



FIG. 1: Density at t = 128 in the simulation with narrow initial spectrum in 960 3 domain. 



seem feasible, and one needs to rely on resolving these 
questions/uncertainties through experiments and simu- 
lations. This article reports a step in this direction. Here 
we raise and give partial answers, based on simulations, 
to the following subset of key questions concerning the 
internal structure of the RTT mixing zone: 

• Analyzing the evolution of h, Rq, and rj with 
time one often observes a non-universal, simu- 
lation/experiment specific behavior, especiall y a t 
transient, so-called early self-similarity, times [la ]. 
Will the relative dependence of scales be a more 
reliable indicator of a universal behavior than the 
time-dependence of the scales? 

• How does the energy containing scale, Rq, compare 
with the width of the mixing layer, hi This ques- 
tion was already addressed in [15|. Here, we will 
elaborate on this point. 

• How different are the turbulent spectra at different 
vertical positions in the mixing zone within a given 
time snapshot? 

• How different are the scales and spectra corre- 
sponding to qualitatively different initial perturba- 
tions? 

The material in this article is organized as follows. We 
start by describing our simulations, we then proceed to 
the definitions and subsequently the observations of the 
various spatial scales characterizing snapshots of the mix- 
ing zone. Finally we discuss self-similarity and universal- 
ity of the emerging spatio-temporal picture of the RTT. 
We conclude by answering the questions posed above. 



II. DESCRIPTION OF SIMULATIONS 

We consider 3D incompressible, miscible Rayleigh- 
Taylor flow in the Boussinesq regime, 

d t v + (vW)v + Vp - vAv = Age, Vv = 0, (1) 
9(C+(«V)c=kAc, (2) 

where k and v are the diffusion and viscosity coefficients, 
while c = (p - PmhO/OWx - Pmin) is the normalized 
density. The Boussinesq approximation for gravity cor- 
responds to fluids with small density contrast, A <C 1, 
where A = {p max - p mln ) / '(p max + p min ) is the Atwood 
number. Here, we restrict ourself to the case of K = v. 

We solve equations ([Hit using the spectral element 
code of Fischer ct al [23| designed specifically for DNS 
of Boussinesq fluids. The equations are solved in the 
nondimcnsional units, 

[l} = (2Ag)-hv 2 3, [t] = (2Ag)-W; 

the results are presented in the same units. This choice 
of units is based solely on the dimensional parameters in 
Eqs. fl][2]) thus reflecting free boundary conditions (ab- 
sence of any wall constraints). The critical wavelength 
and the wavelength of the linearly most unstable mode 
are constants in these units. 

The boundary conditions are periodic in the horizon- 
tal directions and no-slip in vertical direction. The ini- 
tial conditions include a quiescent velocity and a slightly 
perturbed interface between the layers, c(t = 0; z) = 
-6{z + S(x,y)), where 9(z) = | [1 - tanh(0.4z)] is the 
function describing the density profile across the inter- 
face and 5(x,y) is the perturbation. 




FIG. 2: Slices of density and vertical velocity at times t = 32, 64, 96, 128 (left to right) in simulation with narrow initial 
spectrum in 960 3 domain; from top to bottom, the images correspond to vertical slice at y — 480 and horizontal slices at 
z = +0.75/1, 0, -0.75/i. 



We use spectral elements of size 30 3 with 12 collo- 
cation points in each direction. This is equivalent to 
the spectral resolution with the spacing between points 
A = 3. The size of our largest computational domain 
is 1920 x 1920 x 1440 physical units, or 768 x 768 x 576 
collocation points. We stop our simulation at, t = 128, 
when the width of the mixing layer reaches the domain 
size. At the end of simulation, the Reynolds number 
reaches 3? = 7500 to 3? = 13000 depending on the initial 
conditions, where 5i = For comparison, the largest 
Rayleigh- Taylor simulation to date [1 81 ] was performed in 
a 3096 3 domain at resolution A = 1 and reaching time 
t = 248 and 3? = 30000. Our relatively coarse reso- 
lution might raise concerns, especially in diagnostics of 
small structures. Nevertheless, all our results, including 
the spectra and the estimates for microscale rj, are in a 
very good agreement with [l8|, as well as with our finer- 
resolved (but smaller) simulation with A = 1. 

In all cases studied, the (initial) fastest growing mode 
is located at A 24. Most of the presented results were 
obtained in the simulations in the domain of 1920 x 1920 x 
1440 physical units with initial perturbation in the form 
of a narrow initial spectrum, with modes 36 < n < 96 and 
spectral index 0. (Here the spectral index refers to the 
exponent of the wavenumber, as in [24| . and describes the 



shapes of the spectra.) To investigate the influence of ini- 
tial condition we also performed additional simulations in 
smaller domain of size 960 3 (in physical units) with: (1) 
a narrower initial spectrum, with modes 18 < n < 48 and 
spectral index (Figurc[T|); and (2) a broader initial spec- 
trum, with modes 3 < n < 96 and spectral index — 1. The 
two regimes were identified in previous studies as giving 
distinctly different h(t) at transient times (24|. According 
to [U , the first system develops a mode-merging regime 
and exhibits scalings with universal a, while the second 
system develops in the regime of mode-competition, with 
a depending on the amplitude of the initial perturbation. 
We observed that, in spite of the early stage differences, 
the additional simulations gave the same results in the 
turbulent (advanced time) regime as the main set. 



One important focus of our simulations/analysis is on 
resolving the vertical inhomogeneity of the mixing zone 
(Figure \2§ • To achieve this goal we differentiate verti- 
cal slices within a given snapshot, thus calculating vari- 
ous characteristics of the mixing zone such as the energy 
containing scale, the energy spectra and the viscous scale. 
We collect statistics within a given slice z, e.g. contrast- 
ing results for the mixing zone center and its periphery. 
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FIG. 3: Mixing widths h in simulations with narrow (solid 
line) and wide (dashed line) initial spectra. 



the presence of longer wavelengths in the initial spectra, 
and our simulation with broader initial spectrum follows 
the same tendency. As we show below, in spite of the dif- 
ference in a, the two systems are very similar. Further 
broadening of the initial spectrum would require larger a 
and more expensive simulation, while we do not expect 
the results to be significantly different. 

An important thesis of the phenomenology is that 
the internal structure of the mixing zone senses the over- 
all time scale only adiabatically through slowly evolv- 
ing large scale characteristics, of which the mixing zone 
width, h, is the benchmark one. Therefore, our intention 
is to separate the "large scale" question of the overall 
time dependence of the mixing zone width from the set 
of focused "small scale" questions about internal struc- 
ture of the mixing zone. To achieve this goal, we track 
the dependence of the various internal characteristics of 
the mixing zone (see below) on the mixing zone width. 



III. SCALES OF RAYLEIGH- TAYLOR 
TURBULENCE 

A. Mixing zone width 

The mixing zone width is the standard characteristic 
used in the a-studies 0. According to the most recent 
analysis [H|,[IH, the mixing zone width obeys the scaling, 
\fh = -*/ho + t \JuAg. where Kq is an initial-conditions- 
dependent constant. In the simulation with narrow ini- 
tial spectrum, we reproduce this scaling relatively well; 
in the faster-developing simulation with a broad initial 
spectrum the scaling is affected by the finite domain size 
(Figure O. 

The value of a depends on the definition of the mixing 
zone width. Following 0, [ll|, EH we consider the defini- 
tion based on the mixing function, M{c) — 4c(l — c), 

h = J M(c)dz, (3) 

where the overbar denotes averaging over the horizontal 
plane. We prefer integral definitions of the mixing zone 
width over the common definitions based on the values of 
c (for example the half-distance H between two heights 
where c = 0.01 and c = 0.99) simply due to the fact 
that integral quantities are less sensitive to the profile of 
c(z) at the edges of the mixing layer, and consequently 
the quantities are less sensitive to the size of the com- 
putational domain. In the established self-similar regime 
unrestricted by domain boundaries, two definitions of the 
mixing zone width are actually within an O(l) systematic 
factor of each other. The value of the coefficient a, de- 
termined from the slopes of the curves shown in Figure [3] 
at t > 60 are a = 0.029 for the narrow initial spectrum 
and a = 0.040 for the broad initial spectrum. The ob- 
tained values of a are in a good agreement with other 
simulations (see reviews in [9l fl0|). Experimental values 
are higher, a m 0.5-0.7, which is usually attributed to 



B. Energy-containing scale 

The energy-containing scale represents the size of a 
typical turbulent eddy which, intuitively, corresponds to 
the size of the large scale vortices seen in the mixing 
zone snapshot, e.g. shown in Figure Formally, it is 
convenient to define this scale, Ro, as the correlation 
length of the normalized two-point pair correlation func- 
tion of velocity, f(R) = (uj(r)uj(r + R))/(vf), where Vi 
is one spatial component of the velocity vector, v. We 
estimate Rq as a half width of the correlation function, 
/(0)//(i?o) = 2. Defined this way, Rq is consistent (up to 
some 7r-dependent constant) with the wavelength (inverse 
of the wave vector) where the turbulent energy spectra 
achieves its maximum. See e.g. Figure [TT1 

In the developed regime, the correlation length taken 
at the center of the mixing zone grows linearly with the 
mixing layer width, Rq sa h/30 + 7 and Rq ~ ft/17 + 6 
for correlations between horizontal and vertical velocities 
respectively for both types of initial perturbation (Fig- 
ure 5]). In the limit of large i? , this suggests significant 
separation between the two scales, h : Ro = 17 : 1 or 
more. 

Ristorcelli and Clark [l5[ introduced their version of 
the energy containing scale as L = v^/e, and found this 
scale to be of the order of the width of the mixing layer, 
L w OAh. Note that defining L requires a single point 
measurement, while Rq characterizes the two-point cor- 
relations. This single-point nature of L makes it the pre- 
ferred large scale characteristic in the engineering closure 
modeling. Our simulations show that L is significantly 
larger than Rq, where both L and Rq scale linearly with 
h: L f=a 7 — 15i?o ~ 0-5h. The scale separation of L 
and Rq may be viewed as a very favorable fact for the 
engineering modeling of the RTT, thus suggesting a nu- 
merical justification for the closure schemes, e.g. of the 
type discussed in (l~5T |. 
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FIG. 4: The energy containing scale (correlation length) and the viscous scale in the middle of mixing layer in the simulations 
with narrow (solid line) and broad (dashed line) initial spectra. Squares and triangles correspond to scales computed using 
vertical and horizontal components of the velocity, respectively. 
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FIG. 5: Released potential energy and components of the kinetic energy integrated across the mixing layer in simulations with 
narrow (solid line) and broad (dashed line) initial spectra. 



C. Viscous and dissipation scales 



D. Relative dependence of scales 



In our simulations k — v, thus the viscous scale and 
the dissipation scale are equal to each other. (This sim- 
plification reflects our desire to focus on the larger scale 
physics while keeping the resolution domain sufficiently 
large.) We estimate the viscous scale in the middle of 
the mixing layer (z = 0) directly as rj ~ (z^/e) 1 / 4 ~ 
(z/ 2 /(5(Vw) 2 )) 1 / 4 . (The "5"-factor here is an artifact 
of an old tradition, see e.g. [25|.) In magnitude, the 
viscous scale agrees with [18| and with the respective 

phenomenology estimate [HI, r\ ~ ((p/v^h) 1 *. 
The viscous scale decreases slowly with time (see Fig- 
ure|4]). However, our data are too noisy to claim anything 
more than rough consistency with the /i -1 / 8 predicted in 
El El and observed in El El. 



In view of our focus on the internal structure of the 
mixing zone, we choose to study the relative dependence 
of the relevant scales. Thus, Figure [5] shows dependence 
of the energy-containing and viscous scales on the mixing 
zone width. 

Analyzing simulations of RTT corresponding to dif- 
ferent initial perturbations, we confirm the earlier ob- 
served [Til . El ] sensitivity of time-evolution of the mixing 
zone width, scales 77 and Rq, and the integral quantities 
on initial conditions. (See left panels in Figures [Hand [5]) 
However, we also find that the same quantities re-plotted 
as functions of h look very much alike (Figures [4] and 
right panels). Therefore, one conclusion we draw here is 
that the relative scale representation is actually a better 
universal indicator of turbulence within the mixing zone. 
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FIG. 6: Ro - energy containing scale based on the two-point 
correlation function, L - energy containing scale based on 
the single point measurements (see discussion in the text for 
details) , and r\ - the viscous scale, all plotted at time t = 64 in 
simulations with narrow initial spectrum as functions of the 
distance from the middle of the mixing layer. Dashed and 
dotted lines correspond to scale Ro computed using vertical 
and horizontal components of velocity respectively. 



IV. STRUCTURE OF THE MIXING LAYER 
A. Self-similarity 

Figure [7J shows dependencies of the mixing function 
and of the RMS-averaged velocities across the mixing 
layer on the height, z. In magnitude, the velocities are 
0{h}/ 2 ), as suggested by the /i 2 -dcpendence of the total 
kinetic energy: E ~ v 2 h (See Figure[5]). The curves taken 
at three different times are almost indistinguishable from 
each other. This suggests that the mixing zone, viewed 
from the large-scale perspectives, is self-similar. 

Self-similarity of the averaged de nsity and the averaged 
velocities was observed in \m, [H M, M, EE E3 ■ The 
self-similarity itself does not define the specific form of 
the z-averaged profiles, only suggesting that these pro- 
files are smooth functions of z/h. In particular, self- 
similarity is in principle consistent with specific parabolic 
predictions for the size of the dominant eddy and to- 
tal kinetic energy contained in the mixing layer H^j , 
K{z) = K (l - z 2 /h 2 ) and L(z) = L y/l - z 2 /h 2 , where 
Kq and Lq are respective characteristics measured in the 
middle of the mixing layer. Notice, however, that our 
simulation results, shown in Fig. [51 suggest a much flat- 
ter profile for kinetic energy than the parabolic one. 

We notice that when illustrating self-similarity (see for 
instance (l5j). it is common to rescale the mean profiles 
of different quantities using the values at the center of the 
mixing layer and plot these profiles as function of z/h. 
Here, we propose to use /i-based scalings not only for the 
z-coordinate but also for the discussed quantity. Thus, 
in Figure [Jj we rescale velocities with y/h. 

In addition to the mixing function dependence on z, 



Figure [Jj shows 1 - 9{z) = 1 - M(c))/M(c) to be almost 
constant inside of the mixing layer, with a sharp drop- 
off near the edges. Profiles of 9(x), sometimes called the 
molecular mixing fraction, were obtained experimentally 
and numerically in [f| H, 0, [l(| Ht| , and most of the obser- 
vations agree on the fact that at the later stages of the RT 
instability 6{x) remains constant across the mixing layer 
at approximately 0.75-0.8. This suggests that a #-based 
definition of the mixing zone width can be advantageous 
because it generates a more robust scaling. 

The /i-scaling and self-similarity of the mean profiles 
are also observed in the probability distribution functions 
(PDFs) as well as in the correlation functions for density 
and velocities computed at z = (see Figures [5] and O . 
The self-similarity is observed at sufficiently large times, 
t > 64, but is lacking at the earlier times (not shown in 
the Figures). The dynamics one sees at the earlier times 
for the PDF points to transition from a single peak curve 
to two peaks and to a single peak again. The explanation 
for this phenomenon is as follows. The initial, single peak 
distribution is dependent on the initial perturbation of 
the originally sharp interface. The transformation from 
one to two peaks corresponds to transition to the non- 
linear regime of the RT instability, associated with the 
secondary Kclvin-Hclmholtz type shear instability and 
the formation of RT mushrooms. The transition from 
two peaks to one corresponds to the destruction of RT 
mushrooms and formation of the turbulent mixed zone. 
Notice that the emergence of the earlier time transitions 
is consistent with results reported in [f| [l(| Gil for the 
PDF of density, where, probably, the asymptotic self- 
similar regime was not reached yet. 



B. Dependence of scales on z-location 

A discussion of scales is rare in the existing literature, 
not to mention a discussion of how the scales vary across 
the mixing zone. This is partly due to difficulties with 
experimental diagnostics. Notable exceptions arc 15] 
and [ItJ , which discuss the energy containing scale within 
simulations and one-dimensional modeling, respectively. 
The model in (i~7| predicts a parabolic profile for the en- 
ergy containing scale, and this appears consistent with 



earlier time observations in 



and in our setting. At 



later stages the profile in |15[ changes to a much flatter 
one, which also agrees with our estimates of the energy 
containing scale based (as in [l5|]) on single point mea- 
surements. Moreover, we show that later in time these 
profiles stabilize in a nontrivial self-similar solution (Fig- 
ure [10]). 

Our measurements, based on the two-point correlation 
functions of density and velocities as well as estimates 
of the dissipation scale rj, also exhibit monotonic depen- 
dence on z/h. As shown in Figure (6[ both Ro and r\ 
vary insignificantly across the mixing layer with a slight 
increase towards the edges. Among the scales we mea- 
sured, Rq is the only scale which the self-similarity in 
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FIG. 7: Mixing function (left) and velocities (right) vs distance from the center of the mixing layer. The quantities are shown 
at times t — 64 (solid line), t = 96 (dashed line), and t = 128 (dotted line) for narrow initial spectrum. 





FIG. 8: PDF of density (left) and velocities (right) in the center of the mixing layer at times t = 64 (solid line), t — 96 (dashed 
line), and t — 128 (dotted line) for narrow initial spectrum. 
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FIG. 9: Correlation functions of density (left) and velocities (right) in the center of the mixing layer at times t = 64 (solid 
line), t = 96 (dashed line), and t = 128 (dotted line) for narrow initial spectrum. 
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FIG. 10: The energy containing scale based on the single point measurements, L, and based on the two-point correlations of 
density and velocities, -Ro, as function of the distance from the center of the mixing layer at times t = 64 (solid line), t — 96 
(dashed line), and t — 128 (dotted line) for narrow initial spectrum. 



time is still questionable (Figure [TO)) . It might be related 
to a resolution-related systematic error in obtaining the 
correlation functions near zero (Figure [5]). 

Comparing the respective curves for different initial 
spectra, we find that the dependence of Ro(z/h) and 
r)(z/h) on initial perturbation is weak. 



C. Energy spectra 

As expected the level of fluctuations grows with time, 
resulting in a monotonic shift of the turbulence spec- 
tra maxima towards larger wavelengths. The energy- 
containing wavelength, Ao, (corresponding to the max- 
imum) is in agreement with the correlation radius, i?o- 
This applies to various spectra. For example, the maxi- 
mum of the VF 2 -spectrum obtained at z = and t = 64, 
and shown in Figure [TT1 (right), is located at k « 0.15, 
i.e. A ~ 42. The correlation radius at this time is 
i?o ~ 14 ~ X/tt^/2. At time t = 124 the spectrum maxi- 
mum shifts to k = 0.07 (A w 90) while Rq ~ 30. 



Extracting the Kolmogorov scaling for velocity fluc- 
tuations at different wavelengths as well as the largest 
wavelength cutoff (corresponding to the inverse of the 
viscous scale) from the spectral data (e.g. Figure [TT]) 
is problematic due to the lack of spatial resolution. In 
this regard our simulations, as well as many others, e.g. 
@, H, EH EH j lack the power of the simulations performed 
by Cook and Cabot [HI, [l|[ : the latter provide the only 
reliable confirmation (so far) of the expected Kolmogorov 
features of the spectra. (These record simulations are an 
LES run [l]| and a DNS run with 3096 3 points resolu- 
tion [Hj].) Based on our simulation results, wc can only 
state that the range of scales compatible with the Kol- 
mogorov scaling grows with time and that the viscous 
scale decreases with time in accordance with predictions 

of (SGI. 
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FIG. 11: Energy spectra at time h « 180 in horizontal planes z — (circles - horizontal velocities, triangles - vertical velocities) 
and z ± 0.75ft (+ horizontal velocities, x vertical velocities) in simulations with narrow (left) and broad (right) initial spectra. 



V. CONCLUSIONS 

We conclude by presenting a question-and-answer style 
summary for the observations made in the article. 

• Does the relative dependence of the characteristic 
scales constitute a better indicator of self- similarity 
within the RTT than the dependence of the indi- 
vidual scales on time? We found that the energy- 
containing scale, i?o, and the viscous scale, n, ex- 
hibit monotonic evolution with h. At transient 
times both Rq and rj demonstrate much clearer scal- 
ing with h than with the observation time t. 

• How does the energy containing scale, Rq, compare 
with h? We found that at late time, the ratio of 
Rq to h taken at the center of the mixing zone is 
w 1 : 20 and fluctuates little with time. 

• Do the turbulent spectra vary with vertical position 
of the mixing zone ? Analyzing spatial correlations 
for a given time snapshot, we did not observe any 
qualitatively new features of the turbulence spectra 



with transition from the vertically central slice to 
an off-centered one within the mixing zone. 

• How different are the scales and the spectra corre- 
sponding to qualitatively different initial perturba- 
tions? We found that the dependencies of Rq and 
rj, as functions of z/h, on the initial perturbations 
are weak. 

Our effort to extend this type of analysis to account 
for the effects of chemical reactions on the RT Boussinesq 
turbulence will be described in a forthcoming article [28| . 
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fects are in balance, resulting in the so-called Bolgiano- 
Obukhov scaling regime. The 2d phenomen olog ical pre- 
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